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O , The temperature and chemical potential dependent surface tension of bags is introduced into the gas of 

■ quark-gluon bags model. The suggested model is solved analytically. It resolves a long standing problem of 

a unified description of the first and second order phase transition with the cross-over. Such an approach 
is necessary to model the complicated properties of quark-gluon plasma and hadronic matter from the first 
principles of statistical mechanics. In addition to the deconfinement phase transition, we found that at the 
curve of a zero surface tension coefficient there must exist the surface induced phase tranition of the 2"'' or 
higher order, which separates the pure quark gluon plasma (QGP) from the cross-over states. Thus, the present 
model predicts that the critical endpoint of quantum chromodynamics is the tricritical endpoint. 

^ ; PACS: 25.75.-q,25.75.Nq 

^ , Keywords: deconfinement phase transition, tricritical endpoint, surface induced phase transition 

^ 1 Introduction 

The strongly interacting matter properties studied in relativistic nuclear collisions has reached the stage when the 
predictions of the lattice quantum chromodynamics (QCD) can be checked experimentally on the existing data 
and future mesurements at BNL RHIC, CERN SPS, and GSI FAIR. However, a comparison of the theoretical 
results with the experimental data is not straightforward because during the collision process the matter can have 
several phase transformations which are difficult to model. The latter reason stimulated the development of a wide 
range of phenomenological models of the strongly interacting matter equation of state which are used in dynamical 
, simulations. 

ILJ ' One of these models, the gas of bags model (GBM) [3 [HIS]) itself contains two well-known models of deconfined 
and confined phases: the bag model of QGP [5] and the hadron gas model [6]. Hence there were hopes [7] that 
an exact analytical solution of the GBM found in [2] could be helpful in understanding the properties of strongly 
interacting matter. However, this solution does not allow one to introduce the critical end point of the strongly 
interacting matter phase diagram. Also, a complicated construction of the line, along which the phase transition 
order gradually increases, suggested in [7], does look too artificial. Therefore, the present GBM formulation lacks 
an important physical input and is interesting only as a toy example which can be solved analytically. However, 
there are the great demands [8l[9l[T0] for the phenomenilogical models, which can correctly describe the properties 
of the end point of the order deconfinement phase transition (PT) to QGP. 

In statistical mechanics there are several exactly solvable cluster models with the P** order PT which describe 
the critical point properties very well. These models are built on the assumptions that the difference of the bulk 
part (or the volume dependent part) of free energy of two phases disappears at phase equilibrium and that, in 
addition, the difference of the surface part (or the surface tension) of free energy vanishes at the critical point. The 
most famous of them is the Fisher droplet model (FDM) [T^l [131 |T3] which has been successfully used to analyze the 
condensation of a gaseous phase (droplets of all sizes) into a liquid. The FDM has been applied to many different 
systems [H l fTl ] . 

The other well established statistical model, the statistical multifragmentation model (SMM) [121 UHl |T7] , was 
recently solved analytically both for infinite [18l [20j and for finite [211 [22] volumes of the system. In the SMM the 
surface tension temperature dependence differs from that one of the FDM, but it was shown [20] that the value of 
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Fisher exponent tsmm = 1-825 ± 0.025, which contradicts to the FDM value TpjjM « 2.16, but is consistent with 
ISiS CoUaboration data [53] and EOS CoUaboration data 

From the structure of these models, it follows that the GMB can be drastically improved by the inclusion of 
such a vitally important element as the surface tension of the quark-gluon bags. The obtained model is called the 
QGBST model. Its detailed dscussion and the full list of related references can be found in [TTl l25t [26|. 

The great success of the SMM initiated the studies of the surface partitions of large clusters within the Hills and 
Dales Model [27l |28j and led to a discovery of the origin of the temperature independent surface entropy similar to 
the FDM. It was proven that the surface tension coefficient of large clusters consisting of the discrete constituents 
should linearly depend on the temperature of the system [27] and must vanish at the critical endpoint. Thus, the 
Hills and Dales Model [271 [28] is our main guide in formulating the QGBST model. However, for definiteness we 
assume a certain dependence of the surface tension coefficient on temperature and baryonic chemical potential, 
and concentrate on the impact of surface tension of the quark-gluon bags on the properties of the deconfincmcnt 
phase diagram and the QCD critical endpoint. 

Here we will show that the existence of a cross-over at low values of the baryonic chemical potential along with 
the 1** order deconfinement PT at high baryonic chemical potentials leads to the existence of an additional PT of 
the 2"'' or higher order along the curve where the surface tension coefficient vanishes 125". Thus, it turns out that 
the QGBST model predicts the existence of the tricritical rather than critical endpoint. 

The work is organized as follows. Sect. 2 contains the formulation of the QGBST model and analyze all 
possible singularities of its isobaric partition for vanishing baryonic densities. This analysis is generalized to non- 
zero baryonic densities in Sect. 3. Sect. 4 is devoted to the analysis of the surface tension induced PT which exists 
above the deconfinement PT. The conclusions and research perspectives are summarized in Sect. 5. 

2 The Role of Surface Tension 

I begin with the isobaric partition: 

Z{s, T)^JdV exp(-sy) Z{V, T) = prwiW (1) 

L V . /J 

where the function F(s, T) is defined as follows 

oo 

F(s,T) ^ Fh{s,T) + Fq{s,T) = E;=i5.e-'^-0(r,m,) -f u{T) J dv -M-Hs-s^im 

Vo 

At the moment the particular choice of function Fq{s, T) in ([2|) is not important. The key point of my treatment 
is that it should have the form of Eq. ([2]) which has a singularity at s = Sg because for s < sq the integral over 
the bag volume v diverges at its upper limit. As will be shown below the isobaric partition ([T]) has two kind of 
singularities: the simple pole s = s*^ and the essential singularity s — sq The rightmost singularity defines the 
phase in which matter exists, whereas a PT occurs when two singularities coincide [51[TH1[1S]. All singularities are 
defined by the equation 

s* = F(s*,T), (3) 

Note that the exponential in ([2]) is nothing else, but a difference of the bulk free energy of a bag of volume v, 
i.e. —Tsv, which is under external pressure Ts, and the bulk free energy of the same bag filled with QGP, i.e. 
—Tsqv. At phase equilibrium this difference of the bulk free energies vanishes. Despite all positive features, Eq. 
([2|) lacks the surface part of free energy of bags, which will be called a surface energy hereafter. In addition to the 
difference of the bulk free energies the realistic statistical models which demonstrated their validity, the FDM [12] 
and SMM [15j . have the contribution of the surface energy which plays an important role in defining the phase 
diagram structure [18j[22]. Therefore, I modify Eq. ([2]) by introducing the surface energy of the bags in a general 
fashion [20]: 

00 

J V 

Vo 

where the ratio of the temperature dependent surface tension coefficient to T (the reduced surface tension coefficient 

~ T fji"2A;-|-l 

hereafter) which has the form a (T) ^ ^ ■ (fc = 0,1,2,...). Here > can be a smooth function of 
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the temperature, but for simplicity I fix it to be a constant. For fc = the two terms in the surface (free) energy 
of a v-volume bag have a simple interpretation [12j : thus, the surface energy of such a bag is ctqw'*, whereas the 
free energy, which comes from the surface entropy aoT~pV'^, is —TcToT^pV'^ . Note that the surface entropy of a 
w-volume bag counts its degeneracy factor or the number of ways to make such a bag with all possible surfaces. 
This interpretation can be extended to fc > on the basis of the Hills and Dales Model p7t [28] . 

In choosing such a simple surface energy parameterization we follow the original Fisher idea [12] which allows 
one to account for the surface energy by considering some mean bag of volume v and surface v'^. The consideration 
of the general mass-volume-surface bag spectrum is reserved for the future investigation. The power n < 1 which 
describes the bag's effective surface is a constant which, in principle, can differ from the typical FDM and SMM 
value |. This is so because near the deconfinement FT region QGF has low density and, hence, like in the low 
density nuclear matter [35], the non-sperical bags (spaghetti- like or lasagna-like [3S]) can be favorable (see a [5S] 
and references therein). A similar idea of "polymerization" of gluonic quasiparticles was introduced recently [5S] . 

The second essential difference with the FDM and SMM surface tension parameterization is that we do not 
require the vanishing of (t{T) above the CEP. As will be shown later, this is the most important assumption 
which, in contrast to the GBM, allows one to naturally describe the cross-over from hadron gas to QGP. Note that 
negative value of the reduced surface tension coefficient a{T) above the CEP does not mean anything wrong. As 
we discussed above, the surface tension coefficient consists of energy and entropy parts which have opposite signs 
[m [271 [28] . Therefore, cr(T) < does not mean that the surface energy changes the sign, but it rather means that 
the surface entropy, i.e. the logarithm of the degeneracy of bags of a fixed volume, simply exceeds their surface 
energy. In other words, the number of non-spherical bags of a fixed volume becomes so large that the Boltzmann 
exponent, which accounts for the energy "costs" of these bags, cannot suppress them anymore. 

Finally, the third essential difference with the FDM and SMM is that we assume that the surface tension in the 
QGBST model vanishes at some line in fj,B — T plane, i.e. T^ep = TcepifJ-s)- However, in the subsequent sections we 
will consider Tcep — Const for simplicity, and in Sect. V we will discuss the necessary modifications of the model 

with Tcep = TcepifJ-s)- 

The surface energy should, perhaps, be introduced into a discrete part of the mass- volume spectrum Fh, but 
a successful fitting of the particle yield ratios 6J with the experimentally determined hadronic spectrum Fh does 
not indicate such a necessity. 

In principle, besides the bulk and surface parts of free energy, the spectrum ^ could include the curvature 
part as well, which may be important for small hadronic bubbles or for cosmological FT. We stress, however, that 
the curvature term has not been seen in such well established modles like the FDM, the SMM and many other 
systems [13l [14]. A special analysis of the free energy of 2- and 3-dimesional Ising clusters, using the Complement 
method [37], did not find any traces of the curvature term (see a detailed discussion in Ref. [25]). 

According to the general theorem [2] the analysis of FT existence of the GCF is now reduced to the analysis 
of the rightmost singularity of the isobaric partition ([Ij. Depending on the sign of the reduced surface tension 
coefficient, there are three possibilities. 

(I) The first possibility corresponds to a{T) > 0. Its treatment is very similar to the GBM choice ([2]) with r > 2 
[2]. In this case at low temperatures the QGP pressure Tsq{T) is negative and, therefore, the rightmost singularity 
is a simple pole of the isobaric partition s* = sh{T) = F{sh{T),T) > sq{T), which is mainly defined by a discrete 
part of the mass- volume spectrum Fh{s, T). The last inequality provides the convergence of the volume integral in 
([4]) (see the left panel in Fig. 1). On the other hand at very high T the QGP pressure dominates and, hence, the 
rightmost singularity is the essential singularity of the isobaric partition s* — sq(T). The phase transition occurs, 
when the singularities coincide: 



which is nothing else, but the Gibbs criterion. The graphical solution of Eq. ([3]) for all these possibilities is shown 
in Fig. 1. Like in the GBM [21 [7], the necessary condition for the FT existence is the finiteness of Fq{sq{T),T) 
at s = sq{T). It can be shown that the sufficient conditions arc the following inequalities: Fq{sq{T),T) > sq{T) 
for low temperatures and F{sq{T),T) < sq{T) for T ^ oo. These conditions provide that at low T the rightmost 
singularity of the isobaric partition is a simple pole, whereas for hight T the essential singularity sq{T) becomes 
its rightmost one (see Fig. 1 and a detailed analysis of case /^s 7^ 0). 

The FT order can be found from the T-derivatives of sh{T). Thus, differentiating ([3]) one finds 



sh{Tc) 



= Po(Tc) 

m 



(5) 



G + u/Cr-i(A, -cr) • s'. 



Q 



(6) 



1 + uJCr-i{A,-a) 
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Figure 1: Left panel. Graphical solution of Eq. ([3]) which corresponds to a PT. The solution of Eq. ([3]) is shown 
by a filled hexagon. The function F{s,S^) is shown by a solid curve for a few values of the parameter ^. The function 
F(s,<^) diverges for s < sq{^) (shown by dashed lines), but is finite at s = sq(0 (shown by black circle). At low 
values of the parameter ^ = ^yi, which can be either T or /i^, the simple pole sh is the rightmost singularity and it 
corresponds to hadronic phase. For ^ = 3> the rightmost singularity is an essential singularity s = sq(^b), 
which describes QGP. At intermediate value ^ — £,c both singularities coincide sh{£,c) = sq(?c) and this condition 
is a Gibbs criterion. 

Right panel. Graphical solution of Eq. ^ which corresponds to a cross-over. The notations are the same as in 
the left panel. Now the function F(s,^) diverges at s = sq(^) (shown by dashed lines). In this case the simple 
pole Sh is the rightmost singularity for any value of ^. 



where the functions G and /Cr-a(A, —a) are defined as 

OQ 

ICr-aiA,-a) ^ Jdv -p[-Aj-^(^)-^l , (8) 

where A = sh — sq. 

Now it is easy to see that the transition is of the 1** order, i.e. Sq{Tc) > s'^iTc), provided (t{T) > for any r. 
The 2"*^ or higher order phase transition takes place provided s'q{Tc) = s'jj{Tc) at T ~ Tc. The latter condition is 
satisfied when /Ct-i diverges to infinity at T ^ {Tc — 0), i.e. for T approaching Tc from below. Like for the GBM 
choice 1^, such a situation can exist for a{Tc) = and | < r < 2 |^. Studying the higher T-derivatives of sh{T) 
at Tc, one can find a mare general statement, but for our purpose it is not necessary. 

(II) The second possibility, a{T) = 0, described in the preceding paragraph, does not give anything new compared 
to the GBM [2l|7]. If the PT exists, then the graphical picture of singularities is basically similar to the left panel 
of Fig. 1. The only difference is that, depending on the PT order, the derivatives of F{s,T) function with respect 
to s should diverge at s = sq{Tc). 

(III) A principally new possibility exists for T > Tcep, where (t(T) < 0. In this case there exists a cross-over, 
if for T < Tcep the rightmost singularity is sniT), which corresponds to the leftmost curve in the right panel of 
Fig. 1. Under the latter, its existence can be shown as follows. Let us solve the equation for singularities ^ 
graphically (see the right panel of Fig. 1). For (j{T) < the fmiction Fq{s,T) diverges at s = sq{T). On the 

other hand, the partial derivatives ^^"q^''^'^ < and ^^'^q^''^'^ < are always negative. Therefore, the function 
F{s,T) = Fh{s,T) + Fq{s,T) is a monotonically decreasing function of s, which vanishes at s — > oo. Since the 
left hand side of Eq. ([3]) is a monotonically increasing function of s, then there can exist a single intersection s* 
of s and F{s,T) functions. Moreover, for finite sq{T) values this intersection can occur on the right hand side of 
the point s — sq{T), i.e. s* > sq{T) (see the right panel of Fig. 1). Thus, in this case the essential singularity 
s — sq(T) can become the rightmost one for infinite temperature only. In other words, the pressure of the pure 



QGP can be reached at infinite T, whereas for finite T the hadronic mass spectrum gives a non-zero contribution 
into all thermodynamic functions. Note that such a behavior is typical for the lattice QCD data at zero baryonic 
chemical potential [55] , 

It is clear that in terms of the present model a cross-over existence means a fast transition of energy or entropy 
density in a narrow T region from a dominance of the discrete mass- volume spectrum of light hadrons to a dominance 
of the continuous spectrum of heavy QGP bags. This is exactly the case for <j{T) < because in the right vicinity 
of the point s — sq{T) the function F{s,T) decreases very fast and then it gradually decreases as function of 
s-variable. Since, Fq{s,T) changes fast from F{s,T) ~ Fq{s,T) ~ sq(T) to F{s,T) - Fh{s,T) - Sff(T), their 
s-derivatives should change fast as well. Now, recalling that the change from F{s,T) ^ Fq{s,T) behavior to 
F{s,T) ~ Fh{s,T) in s-variable corresponds to the cooling of the system (see the right panel of Fig. 1), we 
conclude that that there exists a narrow region of temperatures, where the T derivative of system pressure, i.e. 
the entropy density, drops down from ^ ^ sq{T) + T^^^^ to ^ ^ sh{T) + T^^^^ very fast compared to 
other regions of T, if system cools. If, however, in the vicinity of T = T^ep — the rightmost singularity is sq(T), 
then for T > T^ep the situation is different and the cross-over does not exist. A detailed analysis of this situation 
is given in Sect. 4. 

Note also that all these nice properties would vanish, if the reduced surface tension coefficient is zero or positive 
above Tcep- This is one of the crucial points of the present model which puts forward certain doubts about the 
vanishing of the reduced surface tension coefficient in the FDM [12] and SMM [15] . These doubts are also supported 
by the first principle results obtained by the Hills and Dales Model [27l [28] , because the surface entropy simply 
counts the degeneracy of a cluster of a fixed volume and it does not physically affect the surface energy of this 
cluster. 

3 Generalization to Non-Zero Baryonic Densities 

The possibilities (I)-(III) discussed in the preceding section remain unchanged for non-zero baryonic numbers. The 
latter should be included into consideration to make our model more realistic. To keep the presentation simple, we 
do not account for strangeness. The inclusion of the baryonic charge of the quark-gluon bags does not change the 
two types of singularities of the isobaric partition llj and the corresponding equation for them ([3]) , but it leads to 
the following modifications of the Fh and Fq functions: 

FH{s,T,^iB) = ^<?je^-'^^^(/.(r,m,), (9) 

OO 

Fq{s,T,^ib) = u(T,hb) dv ^ ; . (10) 



Vo 

Here the baryonic chemical potential is denoted as fiB , the baryonic charge of the j-th hadron in the discrete part 
of the spectrum is bj. The continuous part of the spectrum, Fq can be obtained from some spectrum p(m,v,b) in 
the spirit of Ref. [IHlllS], but this will lead us away from the main subject. 

The QGP pressure pg = Tsq{T, jis) can be also chosen in several ways. Here we use the bag model pressure 

PQ ~ io'^^ ¥ ^ (^)^ 9^^ (^)^ ^ ^^"^ more complicated model pressures, even with the PT of 

other kind like the transition between the color superconducting QGP and the usual QGP, can be, in principle, 
used. 

The sufficient conditions for a PT existence are 



F((sQ (T, = 0) -f 0) , T, = 0) > SQ (T, = 0) , 
F{{sq{T, fiB) + 0),T, ^ib) < sq{T,^b) > ^^A■ 



(11) 
(12) 



The condition ()11|) provides that the simple pole singularity s* ~ sh{T,ijlb = 0) is the rightmost one at vanishing 
/is = and given T, whereas the condition (|12|) ensures that s* = sq{T,^b) is the rightmost singularity of the 
isobaric partition for all values of the baryonic chemical potential above some positive constant ^a- This can be 
seen in Fig. 1 for hb being a variable. Since -F(s, T, hb), where it exists, is a continuous function of its parameters, 
one concludes that, if the conditions Ijlip and ([12)) . are fulfilled, then at some chemical potential ^%{T) the both 
singularities should be equal. Thus, one arrives at the Gibbs criterion ((S]), but for two variables 

SH{T.^l%{T))=SQ{T,^i%{T)). (13) 
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Figure 2: Left panel. A schematic picture of the deconfinement phase transition diagram in the plane of baryonic 
density pB and T for the 2"*^ order PT at the critical endpoint (CEP), i.e. for | < r < 2. For the 3'"'' (or 
higher) order PT the boundary of the mixed and hadronic phases (dashed curve) should have the same slope as 
the boundary of the mixed phase and QGP (solid curve) at the CEP. 

Right panel. Same as in the left panel, but for r > 2. The critical endpoint in the /is — T plane generates 
the critical end line (CELine) in the ps ~ T plane shown by the thick horizontal line. This occurs because of the 
discontinuity of the partial derivatives of sh and sq functions with respect to ps and T. 

It is easy to see that the inequalities ([TT]) and are the sufficient conditions of a PT existence for more complicated 
functional dependencies oi Fh{s,T, pb) and Fq{s,T, ps) than the ones used here. 

For our choice ©, (fTUl) of Fh{s,T, ps) and Fq{s,T, pb) functions the PT exists at T < Tcep, because the 
sufffcient conditions (fTTjl and (fT2)) can be easily fulfilled by a proper choice of the bag constant B and the function 
u{T,pb) > for the interval T < Tup with the constant T„p > Tcep- Clearly, this is the 1"* order PT, since the 
surface tension is finite and it provides the convergence of the integrals ([7]) and ([8|) in the expression ([6]), where the 
usual T-derivatives should be now understood as the partial ones for pB — const. 

Assuming that the conditions (fTTj) and are fulfilled by the correct choice of the model parameters B and 
u{T, Pb) > 0, one can see now that at T = T^ep there exists a PT as well, but its order is defined by the value of r. 
As was discussed in the preceding section for | < t < 2 there exists the 2"'^ order PT. For 1 < r < | there exist 
the PT of higher order, defined by the conditions formulated in [55]. This is a new possibility, which, to our best 
knowledge, does not contradict to any general physical principle (see the left panel in Fig. 2). 

The case r > 2 can be ruled out because there must exist the first order PT for T > Tcep, whereas for T < Tcep 
there exists the cross-over. Thus, the critical endpoint in T ~ pB plane will correspond to the critical interval in 
the temperature-baryonic density plane. Since such a structure of the phase diagram in the variables temperature- 
density has, to our knowledge, never been observed, we conclude that the case t > 2 is unrealistic (see the right 
panel in Fig. 2). Note that a similar phase diagram exists in the FDM with the only difference that the boundary 
of the mixed and liquid phases (the latter in the QGBST model corresponds to QGP) is moved to infinite particle 
density. 

4 Surface Tension Induced Phase Transition 

Using our results for the case (III) of the preceding section, we conclude that above Tcep there is a cross-over, i.e. 
the QGP and hadrons coexist together up to the infinite values of T and/or pB- Now, however, it is necessary to 
answer the question: How can the two different sets of singularities that exist on two sides of the line T — Tcep 
provide the continuity of the solution of Eq. ^7 

It is easy to answer this question for pB < p%{Tcep) because in this case all partial T derivatives of sh(T, pb), 
which is the rightmost singularity, exist and are finite at any point of the line T = Tcep- This can be seen from 
the fact that for the considered region of parameters sh{T,pb) is the rightmost singularity and, consequently, 
sh{T,pb) > sq{T, Pb)- The latter inequality provides the existence and finitcness of the volume integral in 
Fq{s,T, Pb)- In combination with the power T dependence of the reduced surface tension coefficient a(T) the 
same inequality provides the existence and finiteness of all its partial T derivatives of Fq{s,T, pb) regardless to 
the sign of (t{T). Thus, using the Taylor expansion in powers of (T — Tcep) at any point of the interval T = Tcep 
and Pb < I^Bi^cep), one can calculate sh{T, pb) for the values of T > Tcep which are inside the convergency radius 
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of the Taylor expansion. 

The other situation is for /is > ii'^^(Tcep) and T > Tcep, namely in this case above the deconfinement PT there 
must exist a weaker PT induced by the disappearance of the reduced surface tension coefficient. To demonstrate 
this we have solve Eq. ([3]) in the limit, when T approaches the curve T = Tcep from above, i.e. for T Tcep + 0, 
and study the behavior of T derivatives of the solution of Eq. ^ s* for fixed values of /is- For this purpose 
we have to evaluate the integrals /Ct(A,7^) introduced in Eq. Here the notations A = s* — sq{T,^b) and 



= —a{T) > are introduced for convenience. 

To avoid the unpleasant behavior for r < 2 it is convenient to transform 



/C.(A, 7^) ^ g,{Vo) - ^/C._i(A, 7^) + ^/C._«(A, 7^) 



further on by integrating by parts: 

(14) 



where the regular function 5r(Vo) is defined as 



(r-1) 

For r — a > 1 one can change the variable of integration v — > z/A and rewrite /CT--a(A, 7^) as 

exp 



(15) 



VoA 



7 _K 



/C._a 1,7 A- 



(16) 



This result shows that in the limit 7 — > 0, when the rightmost singularity must approach sq{T,^b) from above, 
i.e. A 0+, the function ^ behaves as /Cr-Q(A, 7^) A^"""i + ©(A^""). This is so because for 7 ^ the 
ratio 7^A~'* cannot go to infinity, otherwise the function 1Ct-i (1,7^A~'*), which enters into the right hand side 
of (|14p. would diverge exponentially and this makes impossible an existence of the solution of Eq. ([3]) for T = Tcep- 
The analysis shows that for 7^0 there exist two possibilities: either v = ^'^/S.^'^ Const ot: v = ^'^/S.^'^ 0. 
The most straightforward way to analyze these possibilities for 7 —> is to assume the following behavior 



A = A7" + 0(7^ 



fA^|2[Aa7"-i + 0(7")] 



(2fc+l)Aa7° 
2(T-T„p) ' 



and find out the a value by equating the T derivative of A with the T derivative 
The analysis shows [IS] that for A^^'^ < 77'A^^'' one finds 
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Q-2 



Ai 



an 



2 for T < 1 



Similarly, for A^ ^ > 77'A-'^ " one obtains 7" ^7' 



2fc+l ■ 

' A^^'^ and, consequently, 
for r>l+2^. 



" ~ (r-l)(2fe+l) 

Summarizing our results for 7 — > 0, we can write the expression for the second derivative of A as |25) : 



(17) 



(18) 



(19) 







2fc+l _2 








< 










3-2t 






T-1 
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T < 1 



r > 1 



2/c+l ' 



(20) 



2/0+1 



The last result shows us that, depending on k and fc values, the second derivatives of s* and sq(T,iib) can differ 
from each other for | < t < 2 or can be equal for 1 < r < |. In other words, we found that at the line T — Tcep 
there exists the 2""* order PT for | < t < 2 and the higher order PT for 1 < t < |, which separates the pure QGP 
phase from the region of a cross-over, i.e. the mixed states of hadronic and QGP bags. Since it exists at the line 
of a zero surface tension, this PT will be called the surface induced PT. For instance, from it follows that for 
fc = and K > i there is the 2"'^ order PT, whereas for fc = and k 
order PT, and so on. 

Since the analysis performed in the present section did not include any fj,B derivatives of A, it remains valid 
for the fiB dependence of the reduced surface tension coefficient, i.e. for TcepifJ-s)- Only it is necessary to make 
a few comments on a possible location of the surface tension null line TcepifJ-s)- In principle, such a null line 



i or for /c > and k < 1 there is the 3"^ 
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Figure 3: A schematic picture of the deconfinement phase transition diagram (full curve) in the plane of baryonic 
chemical potential /is and T for the 2"'^ order PT at the tricritical cndpoint (CEP). The model predicts an existence 
of the surface induced PT of the 2"'' or higher order (depending on the model parameters). This PT starts at the 
CEP and goes to higher values of T and/or ^b- Here it is shown by the dashed curve CEP- A, if the phase diagram 
is endless, or by the dashed-dot curve CEP-B, if the phase diagram ends at T = 0, or by the dashed-double-dot 
curve CEP-C, if the phase diagram ends at jis = 0. Below (above) each of A or B curves the reduced surface 
tension coefficient is positive (negative). For the curve C the surface tension coefficient is positive outside of it. 

can be located anywhere, if its location does not contradict to the sufficient conditions pTjl and of the 1"* 
deconfinement PT existence. Thus, the surface tension null line must cross the deconfinement line in the ~ T 
plane at a single point which is the tricritical endpoint (^^^; Tcep{fJ-'^B^)), whereas for /is > /x^^ the null line should 
have higher temperature for the same /is than the deconfinement one, i.e. Tcepi^J■B) > Tc{fiB) (see Fig. 3). Clearly, 
there exist two distinct cases for the surface tension null line: either it is endless, or it ends at zero temperature 
or at other singularity, like the Color-Flavor-Locked phase. From the present lattice QCD data \^8] it follows that 
the case C in Fig. 3 is the least possible. 

To understand the meaning of the surface induced PT it is instructive to quantify the difference between phases 
by looking into the mean size of the bag: 



As was shown in hadronic phase phase A > and, hence, it consists of the bags of finite mean volumes, whereas, 
by construction, the QGP phase is a single infinite bag. For the cross-over states A > and, therefore, they are 
the bags of finite mean volumes, which gradually increase, if the rightmost singularity approaches sq{T, ^b), i-e. 
at very large values T and/or /i^. Such a classification is useful to distinguish QCD phases of present model: it 
shows that hadronic and cross-over states are separated from the QGP phase by the 1** order deconfinement PT 
and by the 2""^ or higher order PT, respectively. 

5 Conclusions and Perspectives 

Here we discussed an analytically solvable statistical model which simultaneously describes the 1** and 2'"* order 
PTs with a cross-over. The approach is general and can be used for more complicated parameterizations of the 
hadronic mass- volume spectrum, if in the vicinity of the deconfinement PT region the discrete and continuous parts 
of this spectrum can be expressed in the form of Eqs. ([9]) and (|10p . respectively. Also the actual parameterization 
of the QGP pressure p = Tsq{T, (ib) was not used so far, which means that our result can be extended to more 
complicated functions, that can contain other phase transformations (chiral PT, or the PT to color superconducting 
phase) provided that the sufficient conditions (fTTjl and (fT2|) for the deconfinement PT existence are satisfied. 

In this model the desired properties of the deconfinement phase diagram are achieved by accounting for the 
temperature dependent surface tension of the quark-gluon bags. As we showed, it is crucial for the cross-over 




dlnF{s,T,iiB) 
d s 



(21) 
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existence that at T — T^ep the reduced surface tension coefficient vanishes and remains negative for temperatures 
above Tcep- Then the deconfinement /j.^ — T phase diagram has the 1"* PT at /is > IJ'Bi'^cep) for | < r < 2 , 
which degenerates into the 2"*^ order PT (or higher order PT for | > r > 1) at fis ~ /^^(^cep), and a cross-over 
for < fiB < IJ-'siTcep)- These two ingredients drasticaUy change the critical properties of the GBM [2j and resolve 
the long standing problem of a unified description of the f** and 2"'' order PTs and a cross-over, which, despite all 
claims, was not resolved in Ref. [7j. In addition, we found that at the null line of the surface tension there must 
exist the surface induced PT of the 2"''^ or higher order, which separates the pure QGP from the mixed states of 
hadrons and QGP bags, that coexist above the cross-over region (see Fig. 3). Thus, the QGBST model predicts 
that the QCD critical endpoint is the tricritical endpoint. It would be interesting to verify this prediction with the 
help of the lattice QCD analysis. For this one will need to study the behavior of the bulk and surface contributions 
to the free energy of the QGP bags and/or the string connecting the static quark-antiquark pair. 

Also in the QGBST model the pressure of the deconfined phase is generated by the infinite bag, whereas the 
discrete part of the mass- volume spectrum plays an auxiliary role even above the cross-over region. Therefore, there 
is no reason to believe that any quantitative changes of the properties of low lying hadronic states generated by the 
surrounding media (like the mass shift of the lo and p mesons [39]) would be the robust signals of the deconfinement 
PT. On the other hand, the QGP bags created in the experiments have finite mass and volume and, hence, the 
strong discontinuities which are typical for the f** order PT should be smeared out which would make them hardly 
distinguishable from the cross-over. Thus, to seriously discuss the signals of the 1** order deconfinement PT and/or 
the tricritical endpoint, one needs to solve the finite volume version of the QGBST model like it was done for the 
SMM [21] and the GBM [22]. This, however, is not sufficient because, in order to make any reliable prediction for 
experiments, the finite volume equation of state must be used in hydrodynamic equations which, unfortunately, 
are not suited for such a purpose. Thus, we are facing a necessity to return to the foundations of heavy ion 
phenomenology and to modify them according to the requirements of the experiments. 

To apply the QGBST model to the experiments it is nesseary to refine it: it seems that for the mixture of 
hadrons and QGP bags above the cross-over line it is necessary to include the relativistic treatment of hard core 
repulsion [40l [41] for lightest hardons and to include into statistical description the medium dependent width of 
resonances and QGP bags, which can, in principal, change our undcrsatnding of the cross-over mechanism [42j . 

Acknov^rledgments. I am cordially thankful to the organizers of the seminar-workshop "New Physics and 
Quantum Chromodynamics at External Conditions" for a warm hospitality and the chance to visit my naitive city 
Dniepropetrovsk and discuss there the physics which is at the frontier line of research. 
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